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iBSTRAOT 


Response of a liglit aircraft with a truss type of 
undercarriage during taxiing and take-off is calculated, 

For taxiing, the excitation is assumed to he stationary 
and the aircraft is idealized to a nultidegree lumped 
parameter model keeping the undercarriage unaltered. During 
take-off, the excitation is treated as nonstationary and 
the aircraft is idealized to a single-degree-of -freedom 
spring-mass system. The response parameters of interest 
are the root mean square value of the displacement and 
acceleration at the cockpit and the dynamic stresses in the 
undercarriage members. The excitation to the system comes 
from the runway roughness. The measurement and power spec- 
tral analysis of runway roughness has also been carried out 
The response of the aircraft during taxiing is computed 
using the power spectral density method and that during 
take-off is computed by numerical integration of the diff- 
erential equation of motion. The effect of tyre stiffness 
on the response is studied. The dynamic stresses in the 
undercarriage members during taxiingare compared with the 
static stresses. The results are summarised in tables. 



NOJ£StTGLATUEE 


JL-J 9 9 • • ’ 

C 

COT 

E 

h 

I 

K, K» 

k 

^ ^2 * * * 

M 


m 

iii-| jin 2? • • • 

P 

P 1 ? P 2 * * ** 
r 

R,(T) 

S 

T 

^ ^2 ^ * * * 


areas of cross-section of undercarriage members 

damping matrix 
covariance 
Young’s modulus 
s amp 1 iiig i nt e rval 

moment of inertia and identity matrix 
stiffness raa,trices 
tyro stiffness 

equivalent stiffness of wing 
lengths of undercarriage members 

mass matrix 

correlation lag value and mass 
elements of M 

displacement vector in norma.! coordinates 
elements of P 
lag value 

autocorrela.tion function) T= rh 
stress matrix 

transformation (direction cosine) matrix 
mode sho,pes 

relative displacement vector in generalized 
coordinates 


X 



X ^ j -Jc 2 ^ ^ 

o 

m 

elementvS of Z 

X 

• 

« 

‘distance along rionway 

y(x) 

m 

« 

ruiway pnoiile 

y(t) 


base acceleration 

z 

* 

• 

absolute displacement vector in generalized 

coordinates 

Z-J » Zg, . . , 

• 

• 

elements of Z 

p 

t» 

• 

crtitical damping factor 

6 

t 

change in length. 

r 

• 

« 

modal participation fa,ctor 

e 

• 

• 

normalized modal cofumn matrix in nomal 

coordinates 

cr 

• 

• 

standard deviation 

2 

• 

mean square relative displacement 

0-2 

• • 

z 

• 

• 

mean square absolute acceleration 

0 ) 

• 

* 

circular frequency (rad, /pec*) 

w. 

# 

• 

"bix 

i na,tural frequency 

n 

• 

• 

spatial frequency (rad,/ft*) 

% «» 

• 

0 

power spectral density of base ‘acceleration 
[ (ft,^/sec.^)/(rad,/ft,) ] 

<tly(W) 

• 

• 

power spectral density of runway rou^ness 
[ ft ,^/(rad,/ft, ) ] 




CONTM^ 

Page 

OEETIEICaTE 



ACKNOWLSEG-EI^ENTS 


abstract 


! 

NOMERCMTURE 



CONTENTS 



CHi-ffTER 1 ; 

INTRODUCTION 

■ 1 

CHAPTER 2 : 

MEilS URISISNT /iND SPECTRilL DENSITY 



CALCULATION OP RUIWaY ROUGHNESS 

5 i 

CHAPTER 3 s 

LUl/IPED PARAIffiTER MODEL OP Ali AIRCRAPT 

21 I 

CHAPTER 4 : 

RESPONSE OP AIRCRAPT TO RTJNY/AI 

' 


ROUGHNESS 

32 

CHAPTER 5 : 

CONCLUSIONS 

59 

REPERET'CES 


62 I 

TABIES 


64 I 

EIGURES 


69 i 



CHjiPTER 1 


CTT RODUOTION 

During taxi, take-off and landing an aircraft is 
subjected to randora dynamic loads which, depend on the nature 
of runway profile (or roughness). The ground induced vibration ; 
environment in which the aircraft operates on the ground is 
significant from the point of view of pilot 's ability to 
read instruments, passenger comfort and structural life of 
landing gear and certain vulnerable locations in the airframe. 
The interest in the behaviour of the aircraft during taxi, 
take-off and landing is thus motivated both by the operational | 
and structura.1 considerations. With the introduction of large, 
flexible aircrafts, the ground loading condition has become i 
more critical. NASA, NACA, RAE and AGARD have undertaken 
several investigations to study the airplane response to runway 
roughness. Roughness of several runways have been measured and | 
: runways have been catagorized as smooth, rough, etc. on rela- [ 
'.ti've basis. The airplane response, such as displacement and | 
acceleration at the centre of gravity of the airplane or at 
the cockpit, to runway roughness during taxiing and take-off 
has been determined mostly through experimental procedures, 

E 

ii 

Hot much work has been done in finding the airplane response 

an „ ■ ; ^ I 

to runway rou^ness through^i.-'ualytical approach, probably 

because of the nonlinear characteristics inherent in the 
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airplane and the randomness of the runway profile. Tung, 

Q 

Penzien and Horenjeff'' hare computed the dynamic response of 
supersonic transport airplanes to runway roughness during 

cl 

taxiing and take-off throu^^ deterministic approach, in which 

the response is computed by numerically integrating the equation! 

of motion of free-free airframe. Nonlinearity of spring and 

damping forces have been taken into account. They have also 

discussed a statistical method to calculate the response of 

airplane to random loading (runway roughness). Recently Kirk 
1 6 

and Perry have carried out the analysis of taxiing induced 
vibrations in a subesasonic aircraft (Boeing 70?) by the power 
spectral density method. The aiipilane has been idealized to a 
three-degree-of-freedom spring-mass-damper system. In their 
work, the excitation coming from the runway roughness has been 
assumed to be stationa,ry and identical on both the main landing 
wheels. The same assumption will be made in the present work 
also . . 

In the present work, the roughness of a runway has I 

been measured and then the dynamic response of a small airplane i 

to runway roughness during taxiing and take-off has been deter- 
an I 

mined through analytical approach. In formulating, the problem 

has been treated as a random vibration problem. The airplane 

under study has a truss t 3 rpe of undercarriage and is idealized 

to a five-degree-of-freedom lumped parsmieter system. The main 
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olDject of tills work is to compute the dynamic stresses in the 
undercarriage memhers developed during taxiing and the absolute 
vertical accelerations at the pilot ‘s portion during taxiing and 
take-off. In computing the response during taxiing, the exci- 
tation to the system is assumed to be stationary. The excita- 
tion to the system comes from the roughness of the runway. 

During take-off at a constant acceleration, the airplane is 
idealized to a single-degree-of-freedom spring-mass system 
subjected to nonstationary input and the response is computed 
by numerical integration of the differential equation of motion* 
The approach to this kind of problem, as will be presented in 
the present work, is new in the sense that random vibration 
approach (power spectral density method) is sought, which seems 
more realistic because most airplanes operate in random vibra- 
tion environment. Also, random vibration studies serve as a 
useful tool in the design and ane, lysis of systems. 

In idealizing the airplane to a discrete or lumped 
parameter model, the undercarriage is kept unaltered and the 
effect of wing flexib:,.Lity is included through an equivalent 
sprung masScAmethod of finding equivalent masses and. lumped 
parameter model of the airplane under study is described in 
Chapter 3. Chapter 2 deals with the measui’eraent and power 
spectral calculations of the runway rou^ness* Two methods 
have been used to measure the runway profile: (i) Surveying 
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, , 1 
and (ii; profilocieter designed by Panchal . The limitations of 

each method are discussed. The runway selected is I.I.T. Eanpur 

airstrip. The .roughness is presented in the form of power 

spectral density functions. Chapter 4 deals with the dynamic 
an 

response of ^idealised aircraft during taxiing and take-off. 

The response is computed in terras of standard deviation or 

mean square value. The conversion of spectral density of runway 

roughness to the spectral density of base acceleration the air- 

to 

plane is subjected ^during taxiing at constant speec’.ls also 
dealt in this chapter. Conclusions are given in Chapter 5, 



CHAPTER 2 


l£BA g.URBMT AIH) SPECTRAL DENSITY 
CALCULATIOH OE RUM AY ROUGHNESS 

2*1, Introduction 

During taxi, take-off and landing the aircraft moves 
over the runway. It has been found that rtinway roughness has 
adverse effects on airplane operation and has posed problems 
to designers. To determine the response of the aircraft during 
taxi, one must have sufficient information about the roughness 
of the runways on which airplanes are going to operate. Such an 
information can be obtained from samples of measured runway 
profile records. This chapter deals with the measurement and a 
brief description of spectral analysis of the runway roughness. 
The roughness is presented in terns of power spectra, 

2 . 2 Measurement of Roughness 

Several methods have been adopted for measurement of 

runway profiles. The conventional and most commonly used method 

is surveying, A brief description of different methods used for 

the runway profile measurement with their limitations is given 
1 1 

elsewhere . Panchal designed and fabricated an equipment, 
called 'Profilometer ‘or 'Profile Eollower’, for measuring road 
surface profiles. His method is faster when compared to other 
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methods and a continuous record of any length can he obtained. 
This equipment was used to get a continuous record of the 
profile of airstrip at I.I.T,, Kanpur in the present work. For 
study of comparison, the runway roughness data was also coll- 
ected by surveying. A brief description of the profile follower 
and its use and limitations in measurement of runway profiles 
ax!e outlined in the following paragraphs. 


2.2.1 Profile Follower 

The profile follower consists of a small mild steel 
wheel of 8 inch diameter which follows the road or runway 
surface profile. The wheel rotates freely on ball bearings 
fixed over an axle the two ends of which are connected to two 
mild steel guide bars, as shown in Fig. 2.1. The guide bars 
slide freely through brass bushes mounted on a channel. This 
channel is used to fix the profilometer to the front bumper 
of a jeep or a motor vehicle. The profilometer is heavy enou^ 
as to nsike contact with the ground when vehicle is in motion. 
The wheel is so designed as to follow completely the profile 
of a minimum wavelength of one foot. 

2.2,2. I ns t rument at i o n 

Main accessories used with the profile follower to 
obtain a record of runway elevation profile were: (i) an 
accelerometer pickup with a vibrationmeter to sense the 
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vibration of tlie profile follower's axle and (ii) a strip 
chart recorder to record the signal transmitted hy the vibration- 
meter. The pickup was mounted on the axle at the place provided 
for this purpose. The pickup at its position senses only the 
vertical acceleration of the axle, the acceleration being 
provided by the runway unevenness. This signal is fed into the 
ibrat ionmet er which integrates the signal twice, giving the 
displacement of the profile follower's axle, which in turn 
represents the elevation of the runway profile, A schematic 
diagram of the instrumentation is shown in Pig. 2,2. An alter- 
nator unit with a voltage stabilizer was used to provide power 
supply to the strip chart recorder. The specifications and 
details of the instruments are given in the appendix of Ref.1, 

2,2,3. Runway Profile Measurement 

The profile follower was fixed to the front bumper 
of a station wagon (jeep) so that it rests vertically on the 
ground. The accessories (vibrationmeter, strip chart recorder, 
voltage stabilizer and alternator unit) were put inside the 
rear portion of the jeep. The pickup was mounted on the axle of 
the profile follower. Electrical connections were made as shown 
in Fig. 2.2, The jeep was driven in the middle of the airstrip 
at a constant speed of 22 m.p.h. Since it was thought that the 
variations in profile across the airstrip were small compared 
to that along the airstrip, only a single record was taken. 
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Tlie irregularities of the airstrip in longitudinal direction 
were of prime importance and hence, as an average, profile of 
the mid-section of the airstrip was measured. The paper speed 
of the recorder was set to 25 mm/sec. 


The criteria for selecting a towing speed of 22 m.p.h. 
of the vehicle are based on two factors: 


(i) the wavelength of interest and 

(ii) the operating frequency (> 2 cps) and 
amplitude ranges of the electronic equipments. 

The variations of runwa,y elevations in wavelengths ranging from 
1 foot to 100 feet v^ere of interest in the present investiga- 
tion, as the rimway profile is made up of superimposed waves of 
different wavelengths, the amplitude generally becoming larger 
as the wavelengths increase. The Y\/avelength of interest, speed 
of the vehicle and frequency of vibration are related by 



( 2 . 1 ) 


where f = frequency (cps) 

V = towing speed of the vehicle (ft ./sec.) 
and X = wavelength (ft.) ■ 

The minimum wavelength of interest (1 ft.) was fixed by the dia' 
meter of the follower wheel. longer wavelengths, according to 
Eq. (2.1), can be achieved by increasing the towing speed of . 
the vehicle or by decreasing the frequency. But, at higher 
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speeds (nore than 25 n.p.h. ) it ?/as found that the vehicle 
vihration would play part in the measurement of runway profile, 
and also the pickup used did not have a flat response at 
frequencies ll^"^ th:'n 2 cps. It worked satisfactorily at 
higher frequencies, from 2 cps to 20,000 cpsi faking all these 
factors into consideration, a towing speed of about 22 m.p.h, 
was chosen, v\7hich would correspond to a maximum wavelength of 
16 ft. 

A sample of the record of runway profile obtained by • 
profilometer measurement is displayed in Fig. 2.3. 

2,2.4. Calibration of Record 

The record of the runway profile on the recorder chart 
paper represents the surface elevations to some vertical scale. 
These elevations are. plotted over a certain mean, this mean 
corresponding to the locus of the centre of the follower wheel 
at every instant of tine. To fix the vertical scale, the pick- 
up was mounted on a shaker table and ?/as connected to the 
strip chart recorder through the vibrationmeter in the labora- 
tory. The connections and settings of switches of instruments 
were consistent with tha,t during profile measurement. The 
shaker table was given a known displacement at a particular 
frequency and a record of the displac ement of shaker table was 
got. The ratio of the displacement on the paper to that of the 
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shaker table was then calculated. This ratio was calculated 
for different operating frequencies of the shaker from 2 cps to 
30 cps. These ratios were found to be constant and an average 
of them would fix the vertical scale. In the present work, 

1 nn elevation or depression on the record represented 1.2 mm 
that on the runway. 

Ik, 

The constant towing speed of the vehicle and the 
recorder paper speed would determine the horizonto,! scale. 
Corresponding to vehicle speed of 22 n.p.h. and paper speed 
of 25 nm/sec., 1 mn on record equalled 1.3 ft. along the 
runway on horizontal scale. 

2.2.5. Limitations on Profile Follower Measurement 

(i) Since the pickup used did not have a flat 

response below 2 cps, variations in longer wavq,'- 
lengths of runway profile could not be obtained 
at low towing speed of the vehicle, 

(ii) At larger speeds (more than 25 m.p.h, ) it was 
found that vehicle vibration was excessive 
causing extraneous vibration of the pickup. 

(iii) Maximum wavelength achievable was only 16 ft, 

though wavelengths upto 100 ft. were of interest, 

(iv) La spite of all efforts to maintain constant 
speed of the vehicle at all instants, there 
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were some fluctuations in the speed and. 
hence lot of concentration from the driver was 
req_uired, 

(v ) Clear records of profiles can he obtained at 
higher paper speed. Unfortunately, hi^er paper 
speed could not be got with the recorder used 
in this work, 

2 *’2 . 6 . Improvement Over Profile follower Measurement 

The assumption of constant towing speed of the vehicle 
could be corrected if the fluctuations in the vehicle speed 
were recorded along ?t/ith the record of the runway profile 
simultaneously as a result of which actual heights could 
be deduced. To achieve this , a 'Make and Break’ contactor 
or 'Timer' (Pig, 2,4) was designed and was fixed to the follo- 
?/er wheel. There was one contact made by the timer for every 
revolution of the wheel. In the present work, the electrical 
signals from the timer and the pickup were simultaneously 
recorded. The timer signals at hi^er towing speeds could not 
be recorded along v/ith the record of the runway profile due to 
the smaller diameter of the follower wheel and the paper speed 
limit. Timer signals could be recorded only at vehicle speeds 
less than 6 m.p.h. With a follower wheel of larger diameter and 
a sufficiently higher recorder speed it would be possible to 
record the timer signals and the runway profile simultaneously. 



12 


2,2.7. Measurement of Runway Profile by Surveying 

To compensate for the limitation (iii) (Section 2.2.5) 
and to compare the results of the measurement made by profile 
follower, the I.I^T,, Kanpur airstrip roughness data was also 
collected from survey measurements. Roughness variations in 
longer wavelengths can be obtained by surveying, Por most air- 
planes the frequency range between G.5 and 35 c.p.s. is the 
^ 10 

region of prime conc ern • At an airplane speed of 50 m.p.h. 
this frequency range corresponds to wavelengths between about 
150 and 2 ft. Since the variations of runwqy profile in 
longer wavelengths of the order of 100 ft. were of interest, a 
reading inteival of 2 ft. was selected which corresponds to a 
minimum wavelength of 4 ft* Variations in runway profile at 
wavelengths less than 4 ft, could not be measured by surveying 
as the staff had a least count of 0,01 ft. Roughness measure- 
ments along the mid-longitudinal section of the runway were 
made by means of a surveyor's level, a tape and a staff. 
Measurements were made for a total length of 1600 ft, of the 
♦ 


runway 
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2.3. Pow er Spectral De a isity Galculation 
Of Rumway Roughness 

The runway profile is a random process with distance 
as the indexing parameter. It provides the basic data for the 
excitation of the aircraft as it moves on the runway. If the 
velocity of the aircraft is loiown, then runway profile gives 
the basic motion as a function of time. Collection of several 
runway profiles thus provides the ensemble for random vibration 
analysis of the aircraft. However, if the runway profile is 
assumed to be ergodic o. single long record can be used. 

The random vibration analysis of the aircraft can be 
performed in t?aD ways; 

1 . Time Domain Approach 

2. Frequency Domain Approach, 

Frequency domain approach is relatively convenient for stationary 
case and is based on the use of power spectral density functions 
(PSDP). PSDP gives the spectral characteristics of the profile. 

As the runway profile exhibits random characteristics which are 
normally ■uniform, the stationarity of the random process is 
assumed in the power spectral calculations. 

Since the runv/ay roughness is a spatial disturbance 
the power spectral density function is defined in terms of the 
frequency argument in radians per foot rather than the 
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conventional argument Oi in radians per second. If y(x) denotes 
the runway height - x is distance along the runway — then the 
power spectral density function, in terms of Q is 

defined by 

(l)j(Q) = lin I / y(x) dx I (2.2) 

where the bars j... | indicate the modulus of complex quantity, 
Eq. (2,2) may be usou to evaluate the power spectral density 
function from observed data, but, in practice, the power 
spectral density may be determined more conveniently and less 
tediously by making use of a related function, the 'auto- 
correlation function’, Ry(T.), defined by 

oo 

Ry(T) = lim f y(x) y(x +T) dx (2.3) 

The autocorrelation function has the symmetrical property 
Ry(7;) = Ry(-T). It is to be noted that the spectral density is 
a non-negative quantity. The autocorrelation and the power 
spectral density functions form a Fourier transform pair and 
are related as follows; 

(fycn) = 

RyC r) = 

—CO 


/ RY(r) e“^^^dT (2.4a) 

mJOO 

OO .Q.t” 

/ (l)y(S2) 


(2A'b) 
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It can -be shown, from Eqs. (2.3) and (2.4h), that 

2 °° 

0 - = y (x) = R^(0) = / (fyCQ) dQ (2.5) 

where cr is the root— necai-sq.uare value (or standard deviation) 
of the disturbance function, y(x), and is a convenient measure 
for comparison of the overall roughness of the runways. 

The evaluation of the power spectra can be done 
through 'data processing' - digital and analog - procedures, 
i Details of measurement and analysis of random data and data 
[processing procedures are given elsewhere”^ In the present 
work, digital technique of data processing was used to evaluate 
the spectral density. The digital technique consists in digit- 
izing a continuous record or data into discrete numbers and 
evaluating the spectral density through numerical calculations 
on digitfil computers. In digitizing, the sampling interval h 
is chosen such that 

b = (2.6) 

where f^ is the cutoff frequency, f^ is also called the 

2 

'Myquist frequency' • The choice of h must be small enou^ so 
2 

that 'aliasing' is not a problem. This aliasing is a source 
of error in which the frequencies higher than 1/2h will get 
folded into the lower frequencies in the range from 0 to 1/2h, 
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The aliasing.' will also confuse the data in the lower frequency 
range. 

In the case of the continuous data of the runway 
irregula,rities obtained from profile follower measurements, the 
runway heights, y(x) , were read off at 0.65 ft. interval, over 
a certain datna, giving a cutoff frequency, f , of 0.7692 cycle 
per f*: or 4.8332 radians per ft. for a runway length of 1600 ft, 
this sampling interval (h = 0.65 ft.) gave a sample size, U, 
of 2462 sample points. In surveying the sampling interval was 
2 ft. giving for leoo-ft. runway a sample size of 800 points; 

The cutoff frequency in this case was 0.25 cycle per ft. or 
1.571 radians per ft. 

The actual steps involved in the numerical calculation? 
of the power spectra are as follov/sj 

2.3.1 . Autocorrelation Eunctions 

For K data values > n = 1, 2, ..., N, the 

autocorrelation function, R^, at the displacement rh is calcu- 
lated from the Bq. (2.3). 

R = R (rh) = 2 y ' y ' r = 0,1 ,2, ... ,m 

r y' ^ R-r ‘'n “^n+r • ’ ’ 

(2.7l> 

where r is the lag number, m is the maximum lag number, and 
y^ is the runway height at a point r, measured over a certain 
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mean. In case of profile follower measurement, was calculated 
from tlie following relation: 


1 




~ ^n ^ ^ 

n=1 


n— 1,2,.,:. ,N (2*8) 


where y^ is the runway height read off over a datum at a point 


n 


, The second term on the ri^t hand side of Eq. (2,8) denotes 


the mean of 




, n = 1 ,2, 


In, the case of survey measurements, it was desirable 
to make some modifications to the actual measured profile, 
because the runway heists exhibited large changes in elevations 
at very long wavelengths. These large changes at low frequen- 
cies have a tendency to complicate and distort power estimates 
at high frequencies because of the effective filter character- 
istics of the numerical estimators. One way of avoiding these 
adverse effects is by *prefiltering' or ‘prewhitening’ the 
original data. The filter used in the present case is a high 
pans digital filter defined by 


y'Cx) = y(x) - y(x - h) (2.9) 

where y‘(x) is the prewhitened runway height to be used in 
Eq. (2.7) for autocorrelation function estimates, y(x) is the 
original runway height measured over a datum and h is the 
internal between successive values of y. The above linear 
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operation (2^9) corresponds to the multiplication in the 
frequency plane hy the function 

F(Q) = 1 - 

where ^ is the spatial frequency in radians per ft. 

The maximum lag mmherj m, is chosen such that it 
is less than one— tenth the sample size N - a thumb rule. In 
the present investigation m chosen was 40, 

2,3,2. Power Spectral Density functions 

(i) The initial or ’raw' estimates of the power 
spectral density function are calculated from the following 
numerical form of Eq. (2.4a) 

'fk = (1) = I [ ^0 + + (4)\ ] 

k = {2*10) 

where (j)j^ = ^ ^ = 0,1,2,...,m. 

^ ^ corresponds to the cutoff frequency, 

(ii) . . Einal ’smooth* estimates of the power 
spectral density are determined as follows: 
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( 2 . 11 ) 

In the case of power spectra calculations for sur- 
vey data, the final smooth power estimates are divided hy 

I PCQ) I - 2(1 - cosQh) 

in order to compensate for the prewhitening operation*. 

Smooth estimates of power spectral density functions 
for "both the data obtained from survey as well as profilometer 
measurements were made with the help of IBM 7044 digital comp- 
uter. Power spectral density fimctions are plotted against 
spatial frequency and are as shown in Pig. 2,5. 

Results 

In Pig. 2.5, the power spectral density function for 
survey data (curve 1) is plotted between frequencies 0,0393 
and 1.571 radians per ft. corresponding to wavelengths between 
160 and 4 ft., and that for profilometer data (curve 2.) is 
plotted between frequencies 0*1208 and 4*8332 radians per ft. 
corresponding to wavelengths between 52 and 1,5- ft, it is seen 
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/froEi "bliese curves that the power estimates from survey data 

are good for longer wavelengths and that from profilometer data 

are good for smaller wavelengths* Though the runway selected 
of 

for two types ^roughness measurement was same, the values of 
the spectral density, estimated separately for each type of 
measurement, seem to agree in the freciuency range of 0,9 and 
1,57 radians per ft. The explanation for the nonagreement of 
values at lower frequencies is obvious from the limitations of 
profile follower measurement, mentioned in Secti on 2«2«5« If 
a slope correction is applied to the data obtained from profile 
follower measurement considering the attitude of the vehicle 
and the follower wheel on the runway, it may be possible then 
that curves 1 and 2 (Pig, 2.5) will agree well in all frequency 
ranges. It can be seen from Pig. 2,5 that curves t and 2 to- 
gether provide a realistic power spectral density curve oyep 

the entire frequency range. This combined KD curve is found 

1 5 

to agree well with the mean PSD curve (curve 5)» 



G^JTER 5 

MODEL OR AN AIRORii FT 

3.1 In trod uction 

The principal structural components of an aircraft 

are : 

(i) lings 

(ii) Ruselage 

(iii) fail surfaces and 

(iv) landing gear* 

The mass and stiffness of the aircraft is distributed over 
each of these components which are elastically interconnected. 

A dynamic analysis of such a, system based on actual mass and 
stiffness distribution is extremely difficult and complicated. 
Therefore, an idealized mathematical model seems necessary to 
dea,l easily with such problems. In this study a typical small 
aircraft v/ith truss-type undercarriage (Eig. 3.1) is idealized 
to a lumped pparrmet or model shown in Rig. 3.2. In arriving at 
the idealized model, the undercarriage is kept as it is, but th( 
half -wings and fuselage (with tail plane) are replaced by 
'equivalent* masses and stiffnesses. 

In this cha.pter, methods of finding the equivalent 
mass and equivalent stiffness for discrete and contlnuotis 
i^systems are described. It is also shown how the idealized model 
shown in Rig, 3.2 is arrived at. 
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5*2. ^Equivalent* Mass and Stiffness 

The motivation behind finding the ^equivalents 
i.iasses ancu stiffnesses (or springs) is the fo.llowings 

It is difficult, sonetimcs, to deal v/ith certa.in 

systems having many or infinite nunber of degrees of freedom, 

Er on finding the equivalent quantities, such as equivalent 

mass, spring, jind applied force, such systems can be reduced 

to systems having one or any convenient number of degrees of 

freedom. The advantage in finding equivalent' quantities is 

that the analysis of the reduced systems can be carried out more 

easily and conveniently by certain known available methods, 

r.nd 

'IDho equivalent mass ^ spring are found by equa,ting the energy 
stored in the equivalent system to the total energy stored in 
the original system, 

5 

V/illiams has given a full description of finding the 
equivalent quantities in his book. The procedure for finding 
equivalent mass and equivalent spring for a discrete system 
having: n degrees of freedom and for a continuous system is 
given below, 'The same procedure ha,s been applied in the present 
study to find the equivalent mass and stiffness of the half- 
wings and the fusela,ge and tail plane.- 

3.2*1. Discrete System 
Example: A Spring-Mass System 

Consider a system having n masses , m^ ,m 2 , . . 



23 


connected by n springs of stiffnesses as shown 

m Fig» 3»3-(a) • Coordinates ^X 2 > represent the displa* 

cefiients of the masses respectively, the notion 

being restricted to the strai^it line through the masses. 
Suppose the equivalent mass, and spring for this system are 
required, which essenti-rally means finding an ‘equivalent^ 
system having one degree of freedom. The simplification in 
reducing degrees of freedom from n to one consists in replacing 
the mass elements of the system by an ‘equivalent na.ss ‘ and the 
springs by an ‘equivalent spring' at the station chosen arbi- 
trarily for the equivalent mass. Once the equivalent system is 
found, standard solution can be ^applied directly. 

Lot the mode shape of the system be defined by ratios 
r. These ratios define the relative displacement of every mass 
in terns of the relative displacement of any two chosen masses. 
Let the motion of relative to m^ be chosen as the basic 
displacement. Then we have 

= r^(x2-x;^), 

x^-x^ = r^(x2-x:^), (5.1) 


If the system is oscillating freely with its centre 
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of gravity at rest and v/itliout any external forces, tlie sum of 
inertia forces is then zero. Therefore, 

+ n2*2 + ^ (5.2) 

Since, in the condition of steady oscillation in simple 
harmonic motion, the accelerations are directly proportional 
to the displacements, Eq_, (-5.2) can he written as 

n^x^ + 112 X 2 + ... + Li^x^ = 0 (3.3) 

Using relations (3.1), Eq. (3.3) can he written as 

n-|X^ + 1112X2 + [ x^ + r2(x2-x^) ] + ... + n^[ x^ + r^(x2-x^) ]= ( 

or 

x^ (n^ + 1:12 + ... + n^) = - (X 2 ‘*x^)(m 2 + + ... + 


or 


x^ = - a(x2-x^) 


(3.4) 


whore 


m 2 + n3n3 + ... + m^r^ 
m-j + m2 + ... +1^ 


(3.5) 


By making use of relations (3.1) and (3.4) the displacements 
of all the masses can he expressed in teicis of the displacement 
of any one mass, such as x^ of m^. Thus* 
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Xg = (t - (s3:j)y 

Xj = (1 - , 

(3.6) 

••«•«• *»••• • 


Y'ihere a is given Toy (3.5). 


Equivalent Mass ; The equivalent mass, M, at station 1 can now 
be readily found by equating the kinetic energy of the equi- 
valent mass to the sum of the kinetic energies of the indi- 
vidual nias'''esi3n this basis, we ha.ve 


1 

2 


H/r* 2 1 / ♦ 2 

MXi == + 


2 


m2X2 




+ m ± ‘ 
n n 


which, by relations (5.6) becomes 


1 2 1 2 2 2 2 \ 

which gives the equivalent mass at station 1 as 

2 2 2 
M = m.^ + m2E.2 + ^®n (3.7) 

Equivalent Spring or Stiffness ; The principle made use of in 
finding the equivalent stiffness is that the same energy is 
stored in the equivalent spring as that stored in all the spring 



of cxie systori for the same refer ere e displacenent , such as • 
Thus, the equiTalert stiffness, K, is given hy 

k2(x2-x^)^ + k^(x^-X2)^ + ... + ^^n^^n""^- 
This , afttvrr using relations (3.6), can ho written as 
•1 Kx^^ = ~ x^^[ k^ + k2(ll2-*'')^ + k3(R5-R2)^ + ... + k^(\-Rn- 
which gives- the equivalent stiffness at station 1 as 
K = k^ + k2(R2-1)^ + k3(R3-R2)^ -f ... + 

(3.8) 

The equivalent system of the present ease is shown 
in Rig;;. 3.3(h) . 

3.2.2. Equivalent Spring-lfoss System for a System 
With Infinite Wupher of Degrees of Ereedom 

Ex.'inple; A Otantilever (Fig. 3. 4) 

In reducing; the cantilever (with infinite number of 
degrrees of freedom) to a single-degree-of-freedom system, we 
have to choose a particular shape for the deflection curve 
(i.e. the node shape), so that the deflection everywhere is 
defined by the deflection at any one point, such as the def- 
lection Jq at the tip, as shown in Fig, 3»4‘* 
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Let the free end of the cantilever be the 'reference 
section' where we v^ant to locate the equivalent mass. If the 
(v3.riable) nass per unit length is n, the kinetic energy of an 
element dx cf the cantilever at x is 

61 = (ni dx) y^ 

1 a 2, ^2 , 

= 2 (^i^) 

where ij = ^ (3.9) 

Yq '^0 

The total kinetic energy is 

T = / 6T = I / n(r^f dx, 

which should be equal to the kinetic energy of the equivalent 
M, Thus, 

whi'ch gives 

1 

M = / in(rt dx (3.10) 

0 

The equivalent spring constant, K, is found by equating- 
the energy stored in the spring under a displacement Yq to 
the total strain energy of the cantilever. Thus, 
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which gives 



v/here E = Young's modulus, and 

I = (Yariable) moment of inertia. 


(3.11) 


It is to be noted that, if the displacement curve 
(mode shra|)e) and the corresponding frequency of vibration are 
determined, then we need calculate only one of the equivalent 
quantities, M ox K, as one can be found out from the other 
uBingthe relation 

c/ - I (3.12) 

where U) is the circular frequency of vibration. 


5.3. Lvtmped Parameter Model of the Aircraft 
Shovm in Eig, 3.1 

In idealizing the aircraft to a lumped parameter 
system, the undercarriage is kept unaltered, but the wing, the 
fuselage (along with tail plane) and the landing wheels are 
replaced by equivalent masses, Ihe half -wings are replaced by 
the equivalent masses m^ and m^ (m^ = i^) and the massless 
members 1 and 12 of equivalent stiffness Kg each, as shown in 
Eig. 3.2. She fuselage with tail plane is replaced by the 
equivalent mass m.^ at the centre of gravity, !I?he main landing 
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wheels are repla.ced hy the equivalent na,sses Ur, and 

^ 4 

(’“^2 ~ ^^4^* tyres are replaced hy springs of stiffness k 
each as shown in the figure. 

Eqs. (3.10) and (3.12) are nade use of in calculating 

the equivalent nass and stiffness for the wing, The equivalent 

quantities are calculated for the fundanental mode, considering 

half wing as a cantilever, The equivalent mass is locaied at & 

station 116 inches from the wing root, l^ith an average weight 

disti’ibution of 12 Ihs./ft, run a,nd average nonent of inertia 
4 

of 60 in. , the equivalent nass, , and the equivalent stiff- 

2 

ness, Kg, for the wing a.re calculated to he 0,4 Ih.-sec, /in. 
j,ind 1319 lbs ./in. respectively. To find the value of equiva- 
lent nasses for the landing wheels, each wheel together with 
undercorrio.go menhers is assumed to wei^ 20 Ihs., giving 
n2 = = 0,052 Ih.-sec. /in. The value of n^ (for fuselage 

and tail plane) is talcen fron the all up wei^t of the aircraft 

(1750 Ihs.) minus the values of , m2, m^, n^. The calculated 

2 

} voliie of rri^ is 3.63 Ihs. -sec, /in. 

The value of the spring constant (tyre stiffness) 
to he used for the landing wheels depends on the tyre pressure. 
Tyre stiffness test at different tyre pressures was conducted 
for an airplane landing wheel having Dunlop tyre - 26”x 7.75* - 
13 " - 8 ply rating. The wheel was loaded vertically in a 
universal testing machine and compression of tyre at different 



30 


loacl.s was noted* Tlie load~def lection curves for the tyre at 
different tyre pressures^ were drawn and are as shown in 
Pig, 3.6, Prom these curves the value of tyre stiffness can he 
obtained. However, these curves could not be made use of in the 
present work as the airplane under consideration had an entirely 
different wheel. Tyre stiffness of 400 , 600 , 800 lbs, /in. will 
be used in the analysis and its effect on the airplane 
response will be studied. 


The idealized model shown in Pig. 3*2 is a symmetric 
structure with masses m^, m 2 , m^, and m^ interconnected by 
members numbered from 1 to 12, The undercarriage forms a 
space structure and consists of members from 2 to 11, Coordi- 
nates X, Y and Z determine the geometry of the structure. The 
generalized coordina,te system for the dynamic model is denoted 
by coordinates z^ , 2 ^, and z^. The following assumptions 

will be made for simplicity in the forthcoming analysis. 


1) Members 1 and 12 are bending members and all 
other members (2 to 11 ) are truss members. 

2 ) All joints, except B and N, are pin- joints, 

3 ) Joints A and M are assumed to lie on the ela- ■ 
Stic axis of the wing so that torsional effects 
are absent. 

4 ) Members 3 to 10 are of hollow round cross- 
sOction and members 2 and 11 are of airfoil 
Gross-section, 
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5) The notion of the aircraft is considered only 

in the -vertical direction (Z-direction) and thus 
the dynanic nodel has only fi-ve degrees of 
freedom. 

The values of various parameters which will be used 
in the response calculations a,re given in Table 3 .I* 



CHAPTER 4 


jjE SPOnSE O P AI RORi iPT TO RUIff/AY RO UG-HHBSS 


4.1. Int r o duct i on 

■During taxiing, take-ofi and landing an aircraft 
is subjected to random excitation due to the runway roughness. 
The iDrccliction of the aircraft response, such as stresses in 
the airfrarie, vertical displacenent and acceleration at the 
c.g. of the aircraft or at the piilot ’s position, etc,, to 
such excitation is of interest. The response can be determined 
either experimentally or anoRytically. In the present investi- 
gation, the response of a light a.irplane to runway roughness 
during*" taxiing, is determined analytically for an idealized 
model (Fig. 3. 2) of the airplane through- random vibration 
approach. In applying this approach, the system is assumed to 
be linear and the excitation stationary. However, v/hen the 
system is nonlinear (nonlinearity comes mainly from landing 
gear system) the analytical approach becomes too difficult. 
Response paraineters of interest are; the stresses developed in 
the undGrearriage members and the vertical acceleration at the 
pilot's portion. These responses are con 5 )uted for different 
constant taxiing speeds and are presented in terms of standard 
deviation or the root-mean-square (r.m.s.) value. Ground 
excitation during take-off is a non— stationary random process. 
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R0spo.tis 0 ^ tak8-*off at "constant acceleration is computed 

Toy numerical integration of tlie equation of motion* In calcu- 
lating the response, the lift forces acting on the airplane 
are not considered, and input to both landing gears are assumed 
to be same. Ihis chapter deals with the analytical procedures 
for finding:* the response of the airplane to random excitation. 

4' • 2 . Response During, laxijng 

An aircraft taxiing at constant speed on a runv^ay 
having a stationary random roughness is subjected to a statio- 
nary excitation. Hence, the laean square response of the air- 
plane can be conveniently determined by frequency domain appro- 
ach which makes use of the |)ower spectral density fimctions 
of the excitation; 

4,2,1, Equations of Motion 

The system considered for response calculations is 
the lumped parameter model shown in Eig. 3*2, Eive degrees of 
freedom considered are defined by the generaliaed coordinates 
and The si»^ of the stiffness matrix for the 
system will be 5 x 5. To get this stiffness matrix, an addi- 
tional degress of freedom at joint &, defined by «here no 
inertia force acts, will be introduced when'll gives rise to a 
6 X 6 stiffness matrix, liim'initing from laae equations of 
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motion will result in 5 x 5 stiffness matrix^ The undemped 
eciuations of motion of the system considering Zg also, can he 
written in the matrix form as 

B‘ir+ K (Z - y) = 0 (4.1) 

where 



(4.2) 


(4.3) 



(4.4) 
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y I0 tlio "ba^ti disp lab 6116111; and ( ) denotes second 'l;iEie*-deri~ 
ya,tive. 

K0W5 to find the total stiffness xnatrixf K, the 
following roiation will he made use of. 

IC = K T (4.5) 

where I is the trfinsfomation riatrix containing the direction 

dosinos of the menihers of the structure referred to coordinates 

rn i 

z^i Z2f T is transposed natrix of T and K* is a 

diagonal natrix containing stiffness of individual memhers. 

Referring to Pig, 5 . 2 , T can he readily written as 


Monbor 

No. 


^2 

Z3 

^4 

"5 

"6 

1 

1 

G 


0 

0 

0 

2 

-h^/L 

2 0 


0 

0 

0 

3 

0 

11^/15 


0 

0 

0 

4 

0 

j 

h^/l^ 

“V ^4 

0 

G 

0 

5 

0 

1 

0 

-N2/I.5 

0 

0 

hg/ljg 

,6 

0 

h^/Le 

0 

0 

0 

f 


7 

0 

0 

0 


0 


8 

0 

0 

CO 

0 

0 

h/Ls 

9 

0 

b 

-h^^/Lg 

^4/^9 

0 

0 

10 

0 

Q 

-114/^io 


0 

0 

i-1 

0 

0 

n-iAi-i 

0 

-Vi'11 ^ 

i 

12 ; 

0 

0 

1 

0 

-1 

G 

13 

0 

-1 

0 

0 

G 

G 

14 ■ ; 

0 

0 

0 


G 

0 
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■vvliLToj aro the lengths of the menhers 

2, 3f ..5 'I’t respectively, h^ = (h^ + h^) and h^ , h^, h^, which 
determine the z-ordina,te of the Tnenhers, are as indicated in 
Pig'# 3 . 2 , It is to he noted that in writing- the direction 
cosines, nemhers 1 and 12 are considered as bending members and 
memberG 13 and I 4 as springs. The stiffness matrix, K*, is 
given by 
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where K,, = equivalent stiffness of half -wing (stiffness of 
riemlDers 1 and 12 ) , 

E = Young’s modulus of elasticity, 
k = tyre stiffness (stiffness of nenhers 13 and 14) 

and ■ii- 2 s are areas of cross-section of the members 

2 , respectively. 

Substituting for K' and T in Eq. (4.8) and 
simplifying the clerdonts.. s'tc., of the total stiff- 

ness matrix arc foimd 'md they are as follows; 

^1 -I = xi 2®^1 /^2 

kgg = 3[ A^/l^^) + Agb^ /Eg ] + ^ 

k^^ = 2Kg + E[ h^ ( 112^^2 ^ "Ml^^ll^^ ^ 

3 3 

+ Ag/I.^) + h^CVl'i + v^4 ■" V^9 ^ 

= E[ + ll^(Ag/I.| + Aic/I'^0^ f 

2 5 

kgg = Kg + A|^b.^/E^-| 

kgg = B[ h|(A5/l| + 1 

ki3 = kj^ = “ 
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^^34 = ^43 = “ Eh|(Ag/Lg + ii^Q/L^Q ) 
kj6 - icg3 = - :sii2(A5/l5 + Aq/L^) 

^46 “ ^64 " E/xya^/L^ 

(4.8) 


and all ivtiior eloniontG are zero. 

. (4.1) is now r-ewrittcn as 

K.. 5 K 


-I'l- 


]m 2 


r°"\ 

\ 


* 

r? -rr 
n 


^21 ; ^22 i / \ 

' J ‘ 


(4-. 9) 


V'/licre 


f \ 



and K^2 » h^^ '^22 partitioned mtrices of K. 

Eqi ( 449 ) is equivalent to the two following equations. 
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- M Z = K^^(Z* - y) + - j) 

The seconu of these equations is solved for zg and 
into the first to pive the reduced form 

- M* r = K*(Z* - y) 


where 


* 


K 






is a 5 X 5 matrix and its elements are given by 


^ i«J ^ ^ 

2 /t 

^22 “ ^22 ^25/^66 

^33 “ ^33 “ ^^36/^^66 
^^44 ~ ^44 "■ ^45/^66 

^55 " ^55 

= k^l = " ■*^11 

^23 “ ^32 ~ ^23 " ^26^36^^66 

^24 " ^42 ~ ” ^26^4(5^^66 

^4 ~ ^43 = %4 ” ^36^46/^66 
* * 

^35 “ ^53 " ^35 ■ 


( 4 , 11 ) 

substituted 


( 4 . 12 ) 
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iill oilier eleEients are ^sero. The equations of motion in the 
ffeneralizGd coordinates are given hj- 


¥r * * ^ At ^ 

M Z + K (Z - y) = 0 


(4.13) 


Putting 


X 


(4.14) 


where 


X = < 


-1 
^2 

X 


4 


X 


5 


(4.15) 


Eq, ( ,-.13) ro.'duces to 


* • • * 

M X + K X 


= - P(t) 


(4.16) 


V'/liere 


P(t) 


in. 






m 




y(t) 


(4.17) 


y(t) is the hase acceleration 
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4 . . 2 . 3^ yncoija.a.od aquations of }Jin±^nr. 

Introc-uoinf prop ort i onal danping into the system, 
the oq'i-'tions notion hecone 

M* X + 0 X + K* X = - a(t) (4.18) 

wl,u.;i'u iJ iii i<hc doniJing nrtrix (5 x 5). These equations are 
st.M couplec;-. , They can be decoupled hy making suitable 

trails L’ernaticms, Let 

# 1 

(4.19) 

Subrititiitin^' this transfornation into Eq, ( 4 . 18) and pre- 

1 

mult ip lying by M ” ? , Eq. ( 4 . 18) bee ones 

_ i 

X’^ + + BX^ - M* ^ E(t) (4.20) 

where r 

A = M*-V2 oM*-V2 

ai'i real) symnetric and positive definite. Ikider these 
conditic'nc and the constraint that /iB = BA, it can be shown 
that thc.ro exists a transfornation 

X = 0 B (4.21) 


Such that 


rn 

G e = I 


(4.22) 
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rcpi’ueunti,! l-ho nomal coordinates, 6 is the normalized modal 
culuriTi I'kitj'ix oj' B, |3 is the critical damping factor assumed 
to be ceiiDt.:.nt in all the nodes and is the i"^^ natural 

.I'rtiCjUi.ncy -f the systen. Trader these transformations Eq,(4,20) 

j’UilUOC.;:'! to 

?!i + + ^.^Pi = ' li ^ = 1,2,.. .,5 

(4.26) 

'.,dn.,re‘ (Y is the i nodal participation factor obtained hy 
sinplifyin{'' the triple matrix product 6 M ' E(t;, The left 
luaid side of Iq. (4.26) can he thought of as representing a 
Binglo-tlerreec-of- freedom systen with natural frequency . 
V/orking back through the transfomations, the mode shapes of 
the system , 112 , ^3’ ^4 ^5’ will he the displa- 

cements in the generalized coordinates, can he obtained. Hence 
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X = U P 

and 

u = H*-V2e 

where 

O' = [ n., tJg Uj U^_ Uj j 

and 


(4.27) 


(4.28) 


“11 


fl'. 

^12 


^15 

f 

I 

'“ 14 " 


■^15 

^21 


^^22 


^23 


"^24 

■ 

^25 

^31 

, Uj = ■ 

j 

^32 

^U3 =< 

U33 

I 

> U, = < 

? 4 

1^34 

'.h=j 

U35 j. 

^41 

1 

I 

^42 

I 


' ""44 

i 

^45 


l"^52j 

I i 

.^^53. 


1 1 
t, 54 J 

1 

.^554 


(4.29) 


are tlie node shapes. 


•i . 2 . 3 . Response G alcu lation s 

If for y(t) a, siraple harnonic forcing function 
y*(t) = e is suDstituted, then the steady-state solution 
of Eq. (4*26) is given hy 


p (-fc) = H, (a>) y(t) 


uv 

1 


(4.30) 


where 


Hp(c4 = 


- 1 




i = i,2, . , .5 

(4.31) 


where again 3‘ = T-1 . Hj_(a!) is the frequency respons 
function of the (linear) system. 


e 
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liie relcitive disp la, cements of the masses in the 
(‘•’eneralj.Jieci coordinates are related to the displacements in 
Ijfiti nor".io,l coordinates by the (4*27) which can be rewri- 
tten as 

5 

>^13 Pj (4.32) 

These displacements are treated as random variables and the 
covariant matrix of the relative displacement vector X, which 
contains the cross products, can be shov^n to be given by^’'^ 


cov 




5 5 

Z Z u 
1=1 j =1 


kl 


u- 




J. 

271~ 


/ I H^Co!) j I 1 dO) 


k a 1 , 2, . . . , 5 
1= 1» 2, 5 


(4.33) 


v/here 

r 2 2 

lipwi = [ 1+ (2?^) ] 


and (OJ) is the power spectral density function of the 
excitation, y(t). 

4.2,4. App roximate Solu tion 

Eor a lightly damped system, jH^(Ct0| given by 
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Eq. (^i-.34) hfeive tegions of pronounoed peaks in ilie hei^Pour~ 
iiood of the oorhesponding natural frequencies If and 
are well separated, the products (Hj_(a)) j }H^ Ik)) jfor i4j are 
seen to he small in comparison v/ith the same products for 
i = Also, terms with 17^3 in Eq, (4 •33) may he negative or 
■positive depending upon the sign of the product Uj^j_ ijfj 

while terms with i=j are always positive. With these approxi- 
mations, Eq. (4.33) may he written as 


00V(3C^X^) = 


f ~ 2 

• IH. ^41 <j>. .(a)) dO) 

2t^‘^ 0 ^ 


k — *1,2, 3' 

1 = 1 , 2,...., 5 


(4.35) 


If the (|Y(Ci)) curve is flat around O)^, then white noise assum- 
ption can he made and the integrals of Eq* (4.35) foi* a- lightly 

'|dcimped system can be approximated by replacing 

Si f 

^discrete values at the natural frequencies Eq.(4.35) 

reduces to 


5 

cov(x|^l) - Ujj-jL u^jj_ ^ 3“ 

■“ i 

k = 1,2, ....5 
1 = 1 , 2, . * • .5 


(4.36) 


In particular, 
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0-2 = .3 

i=i 0i3 Sp 
1 

1,2, ..,,5 (4*37) 

2 

where is the mean square value of the response. 

'1,2,5. Power S pect ral,. Densi ty of Excita tion 

O-'he spectral density, (|)^*-^), of the excitation to 
be used in Eq, (4 .'35) should he in acceleration units like 
■.(f l;.^/!.!0c.^)/(rad./sec.). (|iy (U)) can he obtained from the 
power spectral density functions of the runway roughness in 
'tbfj following way, 

let (^)Y(fi) le the spectral density of the runway 

2 

roughness v.;hich is in ft. /(rad. /ft.) units. An airplane 

■ ■ f, 

taxiing at a constant 'velocity Y(ft./secj) senses an excitation 
of spectral density (|y(<^) ^ ft .V (rad. /sec. ) units, given hy 

( 1 )^( 0 )) = (j)y(Q)/V ( 4 . 38 ) 

Vi/here 

o; = Qy ■ (4.39) 

is in rad. /sec. and Q in rad. /ft. It can he proved from first 



cov(X|^j^) = 


principles that 
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=C0^ (1*y(W) (4,40) 

l.e. [ (j)y(n)/? ] (4.41) 

J^iiua, 1. rora spectral density of the runway roughness spectral 
density of the excitation an airplane senses when taxiing at 
a consl/'i,nt velocity can he found. From the combined spectral 
density curves 1 and 2 for the runway roughness shovi?n in 

2,5 j the power spectral density functions, of the 

base acceleration were computed for different constant 
velocities using Eq. (4,41). The curves of (j)^( 0 ,P) plotted 
against Ci> are as shown in Eig. 4 . 1 , 

4.2.6. '■ndnrf,. Eeviatiq n _or r.m.s. Value of 
Relat ive d isplac ement of Masses 

Q!he covariant matrix, cov(X), of the relative displa- 
ce^mentf.3 in the genera,lized coordinates is a symmetric matrix 
contf;,ining the mean square value of the displacements on its 
diagonal. The square root of the diagonal elements of cov(X) 
gives the standard deviation of the relative displacements 
of masses ,m2, . * & 


4 i 2 i 7 Stan dard Deviation o f Absolute 
Acce leration of fesses 

The absolute acceleration, Z, of different masses 
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''I is given by 

■ 

LO 




i« » • 


Z = X + y(t) 


.'tfi 


V 7 


(4.42) 


Using Bq. (4,27), Eq. ( 4 . 42 ) 


can be written as 


*“ 1/2 .» 

1 e P + y(t) 


( 4 . 43 ) 


Us inf; Eq . (4.22), ( 4 , 43 ) can be written as 


Z 


M* 6 [ P + 0^^ M -/(t) ] 


1 , 0 , 


f f 

Vv'h ere P 


• • 

Z 


U P 


A 




V?.. 


( 4 . 44 ) 


^-P+0 M y(t) gives the absolute accelerations 
in t.h.e i.iorinal coordinates. The mean square value of P^ can be 
fifproxii'iiatcrly ohovm to be given by 


C5' 

Pm 


q .„2 


8 p 


(1 + 4r) i = 1,2, ...,5 


4.45) 


Pa- elements of P^. Neglecting the cross 

Xn'oducts when the frequencies are v/ell separated and using the 
orthogonality condition in the normal coordinates that the 
covariance of p. p. is zero, the mean square value of the 
absolute acceleration of the masses m^ , 1 ^ 2 , m^ is given by 


• *2 

z. = 

^ D=1 


2 2 .*2 

s U.J PA. 


X — ••#5 5 


(4.46) 



the sqUcire root of which gives the standai-d deviation of the 
-■■hsolute acceleration. Incidentally, the absolute acceleration 
of the mass (I'ig. 5 .2) represents the absolute acceleration 
at the cockpit. 


4- . 2 , 8 . Stand ard Devi atio n of Stress in 
Underc arri age Members 

fhe dynamic stresses induced in the members 2 to 11 
(soG Pig., 3 .2) during taxiing are of prime interest here. Let 
"^11 ^ud s 2 } j ^*^1 be rhe changes xn 

len^h and the stresses developed, respectively, in the 
members 2,3»,.*»11 and be denoted by 



(4.47) 


and 



(4.48) 
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6 ir. rulciotd bo the relative dieplac einent vector X by 

6 = 1 X (4449) 

where i is the transfonnation matrix given by after 

deleting the first and the last three rows. Ogj 6 ^, ...5 
.•.'iT'o dependent on each other. The covariant matrix of 6 is 
given b>' 

cov(6) = T cov(X) ( 4 . 50 ) 


cov(X) in Eq. (4.50) is a 6x6 symmetric matrix and includes 
the relative displacement, Xg, at the coordinate Zg. Xg can 
be obtained from the last equation of the equations of motion 
given by Eq. (4.1) as 

Xg = Zg-y = - (kg^x^ + kg2X2 + kg3X^ + ^64^4 ^65"'5^/^66 

(4.51) 

The cov(XgX^), i = 1,2,. ..,6 can be readily found by using 
Eq. (4.51) and the cov(x^x^), i, 3 = 1 ,2, . . . ,5. Now, the 
changes in length given by 6 are related to the stress matrix; 

by 

S = K' 6 ( 4 . 52 ) 

where K' is a diagonal matrix given by 
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E' = rVl-iJ (4.53) 

•Ihc! c ovtU'iciij, b Liatpix of S is givoii "by 

cov(S) = K' cov( 6 ) ( 4 , 54 ) 

•Th© fc.(,j u..,.u'fci root of tlie diagonal elements of cov(S) gives the 
Eitcindcird devio-tion of the dynamic stresses* These, stresses 
denott: the standard deviation of the fluctuation of stresses 
nbc^ut idle static stresses and may be added to the static 
,'i'tresuc's to get estimates of meiximum dynamic stresses, Eor 
oxanip.Lo , stress given by j^ug [ + 30 p represents an estimate 
of mc-iximum stress with exceedance probability less than 0.27 
perotait;, vs/here is the static stress and 0 ^ is the standard 
deviation of the dynamic stress, 

4,2.9. Result s^ 

The natrual frequencies and the mode shapes for 
the idealised system of 'Piper' airplane have been computed 
using a standard subroutine in IBM 7044 computer. The values 
are as shovm in Table 4 .I. The standard deviation of the 
absolute acceleration at the cockpit, and the dynamic stresses 
in the undercarriage members for different taxiing speeds 
have been computed and the values are as shown in Table 4.2. 
The probabilities that the dynamic stress will exceed the 
root mean square value o^, 2a^ and 3a^ are given in Table 4.3 

' llT.KAHPm. I 


I 



CT 
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Thu ivitios of the estimated maxiraum dynamic stress to the 
St otic stress have been calculated and are shovm in Table 4.4. 

ihe effect of tyre stiffness on the fundamental 
frequency and the response calculations for the idealized model 
haS) been s ladled. The fundamental frequency and the response 

.tound 'bo increase with the increased value of tyre stiffness 
The vai’iation of the fundamental frequency and the absolute 
accul oration at the cockpit against tyre stiffness has been 
shown in Table 4.5. 

Rospo n se During Take-off 

Many practical problems in structural d 3 ''namics 
involvti excitations which are nonstationary in character. For 
exoinplCj an airplane experiences nonstationary excitations 
during take-off end landing. These nonstationary excitations 
l€>ad to nonstationary responses. In practice, the analysis 
of such problems involving nonstationary excitations is diffi- 
(Milt to deal with. Also, such problems require many field 
records for estimating the statistical properties of non- 
stationary random processes. In the present study, the response 
of an airplane (Piper) to nonstationai^,^ input during take-off 
is determined through a deterministic approach, in which the 
response is computed by numerical integration of the differ- 
ential equation of motion. The airplane is assumed to move at 
constant acceleration 'a’ during taJee-off. 
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4 • y] • solution of tlis P rnhi ptt^ 

Iho a.irplane under study is idealized to a single-^ 

dogree-of “freedom system shown in Pig. 4*2, The mass m is so 

chosen as to represent the all up v/ eight of the airplane 

(1750 Ih.), The system is assumed to have a natiiral frequency, 

(jt^f equal ‘bo the fundamental frequency calculated for the 

lumped parameter sj'-stem (13*27 rad. /sec.) and a stiffness 
2 

lo szO) m. Introducing viscous damping c, the equation of motion 
for thii'j idealized system subjected to a base displacement, 
y(t), is written as 

m ‘3 + c (z - y) + k (z - y) = 0 (4*55 ) 

where (* ) denotes differential v/ith respect to time. Dividing 
by m and setting c = 2p6)m and Ct)= l/k/m, Eq, (4*55) canbe 
written as 

+ 2(30) z + 0 / z = 2pCDt +c/y (4.56) 

where p is the critical damping factor and co is the 
natuial frequency of vibration of the system* 

Assuming that y may be approximated by a segmentally 
linear function as shown in Pig, 4.3, Eq, (4*56) may be 
written as 

P 2 r -^i / ^ \ 1 

V+ 2f3ai^ +60 Z = 2pa>75^ + a> [ Ji + 5E7 j 

< t < 


( 4 . 57 ) 
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v/liere 


and 


Aji = 



(4.58) 


Thu solution of (4.57) ? ^ t ^ "^ 1+1 glTen by 


z = ’[ sin6t)l/S-p^(t-t^) + cosa>/l-^ (t-tf) ] 

+71+ ^ (t-ti) (4.59) 

in which and C 2 are the constants of integration. Setting 
z = Zi and S = Zi at t = t^ and solving for and O 2 J it 
can bO’ shown that 


= (i/tdi/i-p^) [ Zi + P<^.4^(xi - yi) -Aji/At^ ] 

( 4 . 60 ) 

and C 2 = ^^1 “ y^. 

Substituting (4.60) for and C 2 in Eq. (4.59) > z and ^ 

at t = t. j can be shown to be given by 
i+1 


z^^^ = A(p,C4 At^) Zi + B(pM A"tl) ^1 


(4.61a) 


in which 
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^i+r 


^i+1 


, z . 

’ 1 


Lz 


i+1' 


2 * 




( 1 


3 y*! *“ 


14- 1 


(4.61b) 


A 


^"11 ^"^ 12 " 



'bi 

^12 


J 

B = 



_ ‘"^21 ^22 _ 



f 21 

^22 _ 


(4.61c) 


The elements of matrices A and B are given by 


‘"'ll 


e sin (j) i/l-p^ At. + cosCi>l/l«P^ At.) 

fi-r ' ^ ^ 


a 


pa>At . I — -2 

^o - ■ sinOJtl-p At. 

I c. to ' ‘ 

aivi-^"^ 


(4.61d) 


a 


21 




^ sina;i/i-^^ At, 


cl 


22 


^-pwzsfc^ (cosWlf^^ At. - -p^o sinOA^I^ At . ) 

^ fi-r 
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h" " 

c^Vi-P At, 


(4.61e) 


r _ 1 • 


i [ ~~=^ (O)- l:p-) sinCO'/T^^ At . + 


ffi- <=oswVl-P%^t. ] - 1 


1 


bgg = o-l*^l 


[ •p sJ. w— - S jj3ft>')/l-p^ At . 

■/rip'^At. ^ 


— ooeCwfl^^ At 
. 1 


Prom Eg. (4.56), the absolute acceleration, *z' 5 _, of the, mass 


m o.t time t^ is fiven by 


= 2^.Qj(S-j_ - Sj_) +w^(yi - Zi) 


(4.62} 


(Zi - y^) and (z - y) in Eg. (4.62) respectively denote the 
relative velocity and rela.tive displacement of the mass m. 

Hence, if the displacement and the velocity of the system are 
known at some time t^, the state of the system at all subsequent 
times t^ can be computed exactly by step-by-step application 
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oi L(],,, . (4.61) .and (4.62). The computational advantage of this 

approach lies in the fact that the matrices A and B depend only 

on p, (j) andAt^. If At^j the interva.1 of digitization, is 

constant (that is, the runway profile record is digitized at 

equal intervals), the computational efforts will still be 

sirapler in the sense that the matrices A and B which consist of 

big expressions [ Eqs, (4.61d) and (4.61 e)] need to be evaluated 

8 

only once at the beginning of the response calculation In 
the pros ant investigation, the response is computed for the base' 
displacement y^ measured at equal time intervals, Ab^* The 
interv;,'! of integration, AT, is chosen to be one-fifth of the 
interval of digitization At^..- The time interval of digitization 
is obtained as follovi/s.- 

If the airplone, initially at rest,- is moving at a 
constant accolerction 'a') then the total time tahen to cover 
a distanc e ’ s ' is given ' by 

t ’ ■ = '/Bs/a 

The time interval of digitization is given by 
t = tyif 

Where N is the total number of sample points. 
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4 . 5 *?. 

'i-'lit; :t‘( oponse of the airplane moving with a constant 
aceoleration 'a’ during take-off is computed for a take-off 
distance of s = 300 ft. The runway profile is digitized at 
2 fl;. interval which gives N = 150. The response quantities 
c'; 0 'n(tut(;il 'U'c; tlio maximum and the root-mean-squa.re values of the 
.nolativi' displacement and the absolute acceleration of the 
muss m, which in a way respresent those at the cock-pit. The 
rc.'sult!..! ;iiu' suirmarised as shown in Table 4.6. 



CHaPT BR 5 

GONamsioN s 

The thesis consists of two parts: 

I. Measurement and calculations of povv?er spectral 
density of runway roughness, 

11, Responf^e calculations of a light a,irplane to 
runwa,y roughness during taociing and to,ke-off. 

The measurement of runv/ay roughness (l.I.T., Kanpur airstrip) 
has "been carried out by two methods: (i) surveying and 
(ii) profilometer designed "by Panchal . The power spectral 
calculations of the roughness have been carried out on a 
digital computer, from the power spectral density curves 
(Rig. 2.5) it is observed that the roughness data obtained 
from surveying is good for longer wavelengths (4 ft. to 160 ft.) 
v/hereas the data from the profilometer measurement is good 

for smaller wavelengths (1.3 ft. to 16 ft.). The limitations 

♦ 

on the range of wavelengths mea,surable by the profilometer are 
dictated by the characteristics of the pickup and the towing 
speed. It is foimd that at speeds greater than 25 m.p.h. the 
vehicle vibration becomes excessive. It is suggested that a 
multiwheol mounting for the profilometer can be used to 
measure the roughness in longer wavelengths. The power spectral 
density curve obtained for the airstrip is found to agree well 

15 

with the HATO curves. 
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-Pile response of a light aircraft (Piper) Wi1;h a truss 
type 01 undercarriage during taxiing has been computed by the 
power spectral density method, lEhe a,ircra-ft ha .q been idealized 
to a f iTO-degree-of-freedom lumped parameter system. The 
response panrmeters computed are the r.m.s. va-lues of the 
absolute acceleration at the cockpit and the dynamic stresses 
in the undercarriage members. The prima-ry contribution to the 
aiipl'-aie rosponse is found to be from the fundamental mode only. 
The response is found to increase as the taxiing speed is 
increased. The estimated maximum stresses in the members, vi/ith 
an exceedance probability of 0.27 percent, are found to be 
betv/een 1.48 and 1,93 times the static stresses. The tyre 
stiffness is found to have a significant influence on the 
frequency and response calcula.tions . The response increases as 
the tyre stiffness is increased. The fundamental frequency of 
the aircraft under study is found to lie betv/een 2,112 and 
2.983 cps for a tyre stiffness between. 400 and 800 lbs. /in. 

The absolute acceleration at the cockpit is found to lie 
betwoon 0,150 g and 0.277 g for a tyre stiffness between 
400 and 800 lbs. /in. at a taxiing speed of 20 m.p.h. 

Por the response during: take-off the airplane has 
been idealized to a single-degree-of-freedom spring-mass 
systa:;!. The excitation in this case is a non— stationary random 
process. The response (absolute acceleration at the cockpit) 
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is conijiit ed by Hums ric ally integrating the eq.uatioii of notion, 

For B. tci-ke— off speed of 40 m.p.h* and a take— off distance of 
300 ft,, the r»m,s, vri.li,ie a.nd the maximum value of the absolute 
acceleration at the cockpit are found to be, respectively, 

0.2^1 g and 0.369 g. 

It can be concluded that the acceleration levels 
computed for the small aircraft during taxiing and take-off 
aro Vtfitliin the tolerable limits. But, for 1^3?ger flexible air- 
crafts the response will be critical to runway rou^mess 
during tcixiing and take-off* The technique presented in this 
work can be used for determining the response of large air- 
crafts to runwo.y roughness , during taxiing and takefr»off. It will 
be interesting to calculate the response of aircrafts to 
runway roughness during landing also. The response during 
landing can be ca,lculated by the same method used for calcu- 
lating the response during take-off, i,e, by the numerical 
integration of the equation of motion. On the basis of the 
calculations of dynamic stresses presented in this thesis, 
the fatigue life of the undercarriage structure can be estimated. 
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TiiBIE 4.1 ; EaTUSAI P i{EQlJMGrES MliD MODE SHAPES OP SHE 
idealized model (EIG. 3.2). 

(CASE; LYEE Si'IPEHESS k = 400 lbs, /in.) 
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TABIE 4.3. s PROBIBILHY OP iSCCEEIjiNOE OP DYElI^flC STRESS 


Stress 

°I) 

pPo|, 


j Prob. of 

1 Exceedance i 

'31.73 

I 4.55 

0.27-/ 

* ♦ 


* 

TASTE 4.4- FulTIOS OP ESTnETTED IvlAXETUIT PYl'UMIC STRESS TO 

STATIC STRESS (CASBs TYRE STIPPNESS, k = 4001Es^/in.) 


Member lo. 

2 I 3 

I 

4 1 5 

; 1 

f 

6 I 

i 

7 i 

8 

9 1 10 

\ 

j 

11 

t 

Ra.tio; I 

I 



1 


i 



1.55:1.38 

i i 

1.38, 1.68 

i ‘ 

^ 1 

1.38 

r ,j 

1 .38 

[ I 

1.68 

L : 

1.38|1.38 

j 

1.55 

* 20 

j ; 

1 .70i 1 .48^ 

1 .48i 1.93 

il.48 

il.48 

1 . 

|1.93- 

1.4811,48! 

:1.70 

i : 

; j 

i ) 

! i 

! ' ’ 

} i 

i ' 


! 

i 

1 

i i 

i ! 

! 1 ^ 



* Estimatea Ma:-c. Dyn. Stress = j pg ]• + 3a^ 



TABLE : 


4.^:. - EPTEOT OB TYRS STIFBlirESS ON BU..)]DAI£ENTAL 

brequm'Icy ailj 'absolute acceleration at the 

COCICPIT (CASE; TAXIE'G SPEED = 20m.p.li.) 


Tyre stiffness 
lbs ./iii. 


400 

SOO 

800 


j ~ . 

I Eundamental 

Abs. Accln. 

I Frequency (cps) 

at cockpit 

J j 

- 

2.112 

i. 

0.150 g 

I 

2.585 ! 

i 

0.213 g 

I 2.983 

0;277 g 



TABLE 4, '6 : RESPONSE DURB'TG TAiCE-OPE 

I ',>!**’ ■ ■■ 

( 1 . Rsla;bd,T?'e Displa-cement of mass rn: 

(i) r»m,s. Ta,lue = 0.3425 in. 

(ii) msiX. value = 0,8228 in. 

1 

2. Absolute Acceleration of mass m; 
(i) r.m.s* value = 0.221 g 

I (ii) max* value = 0.369 g. 
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FtG2.5_POWEH SPECThAL DENSITY OF RUNWAY 
illlTilliiiTllE 








(b) 

FIG, 3.3(a) SPRING -MASS SYSTEM ( n d.o.f.) 
(b) EQUIVALENT SYSTEM ( i d.o.f ,) 





WING .. ("0= 












FT. /5EC.‘*)/(RA; /SEC 





■' :*^'^Frequency.,,w(RGd;^5<3c.; 

FiG',_4.1 ^ POWER S'pECTRAC ? 'DENSITY OF 
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FIG. 4 2 _ SINGLE - DEGREE -OF - FREEDOM 
system: 


^.pSpur, 


FIG.4..3_ IDEALIZED BASE DISPLACEMENT 



